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ABS'l’RACi’ 


Total  electron  content  data  (TEC)  have  been  deduced, 
by  the  Faraday  rotation  method,  from  the  signals  of  the 
Spanish  beacon  satellite  lilTAoAT  for  the  period  Feb,  1975' 
Sept. 1976,  The  starting  date  of  the  period  was  establis- 
hed after  an  analysis  to  determine  the  time  when  the  ef- 
fect of  the  resid\ial  spinning  of  the  satellite  was  negll 
gible  as  compared  with  the  Faraday  effect.  Monthly  mean 
TEC  values  have  been  obtained  from  these  data  for  evety 
latitudinal  degree  of  their  sublonospherlc  points  for 
six  half  an  hour  intervals  centered  at  0900,  0930  and 
lOOOLT  in  the  mpming  and  1930,  2000  and  2030LT  in  the 
evening.  These  monthly  values  have  been  analyzed  to  ob- 
tain seasonal  and  latitudinal  TEC  variations,  and  a mo- 
del has  been  deduced,  that  gives  TEC  values  as  a func- 
tion of  the  latitude  and  the  month  of  the  year  for  the 
range  of  latitude  of  the  Iberian  peninsula  and  the  intSr 
vals  stated  above. 
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IKTRQDUGTION 

Polarization  angles  of  the  signal  of  the  Spanish  bea- 
con satellite  INTASAT  have  been  recorded  at  the  Obsserva-  I 

torlo  del  Ebro  for  the  whole  active  life  of  the  satellite » 
from  Nov,  1974  to  Oct,  1976, 

Due  to  the  fact  that  the  satellite's  orbit  was  Sun-sin 
cronous,  the  records  correspond  to  a short  period  In  the 
morning  (between  0830LI  and  llOOLT)  and  another  short  pe- 
riod In  the  evening  (between  1900LT  and  2100LT),  Before  re 
duclng  the  data,  an  analysis  has  been  done  to  determine 
the  time  when  the  satellite  was  properly  stabilized  In  or- 
der to  obtain  reliable  data  of  total  electron  content(TEC) 
by  the  Faraday  method. 

About  1600  passages  of  the  INTASAT  have  been  analyzed 
to  obtain  seasonal  and  latitudinal  TEC  variations  and  a nu 


merlcal  model  has  been  deduced,  that  gives  TEC  values  as  a 
function  of  the  latitude  and  the  month  of  the  year  for  the 
time  interval  covered  by  the  data, 

SATELLITE  STABILIZATION. 


As  Is  well  known,  the  accurate  obtentlon  of  TEC  by  the 
Faraday  method  requires  a good  stabilization  of  the  satelli 
i te  that  eliminates  its  spinning  about  an  axis  different 

from  tha  transmitting  dipole,  V/hen  this  condition  is  ilot  ao 
' compllshed,  the  Faraday  effect  appears  In  the  records  mixed 


up  with  variations  due  to  the  orientation  changes  of  the 
transmitter  antenna. 

In  the  absence  of  Faraday  rotation,  the  signal  transmit 
ted  by  a dipole  that  rotates  about  Its  center  generating  a 
cone,  Is  received  in  a fix  antenna  with  the  strength  varia- 
tions shown  in  figs,  la-lc,  depending  on  the  relative  posl- 


a 


b 


c 


Fig.  1,-  Strength  signal  variation  produced  by 
a rotating  dipole  on  a fixed  antenna 
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tlon  of  the  two  antennas  (of.  Roger  and  Thomson  1960) .When 
the  receiving  antenna  is  parallel  to  a line  inside  the  co- 
ne generated  by  the  transmitter,  the  received  signal  corres 
ponds  to  fig.  1%  if  the  receiving  dipole  is  also  Inside  the 
cone  and  to  lb  if  it  is  outside.  Fig.  Ic  corresponds  to  the 
case  when  the  dipole  is  perpendicular  to  the  previous  cases. 

The  effect  of  the  interaction  of  satbllite  spin  and  Para 
day  rotation  is  two-foldt  a)  the  number  of  fadings  increa- 
ses, and  b)  the  variations  on  the  receiving  antenna  change 
from  fig.  la-b  to  fig.  Ic  and  the  other  way  roiuid,  because 
the  rotation  of  the  polarization  plane  due  to  the  TEG  is  equ^ 
valent  to  rotate  the  transmitting  dipole.  This  effect  is 
clearly  seen  in  fig.  2 v/here  the  record  corresponding  to  the 
orbit  120  of  the  li'fTAOAT  is  shown. 

The  effect  of  the  anomalous  spin  was  gradually  diminis- 
hing for  several  days  and  finally  disappeared. 

Several  authors  (Aitchison  et  al.  1959)  discussed  the  ef 
feet  of  the  interaction  of  Faraday  rotation  and  satellite 
spin  to  deduce  aproximate  values  of  TEC  from  records  of  two 
separate  frequencies  (90  ISHz  and  40  MHz).  We  intended  to  de- 
termine a date  from  which  the  effect  of  the  satellite  spin 
was  negligible  as  compared  with  the  Faraday  rotation  in  or- 
der to  obtain  acctirate  values  of  TEC.  With  this  pupose  we  ma 
de  a Fourier  analysis  of  the  signal  amplitude,  recorded  for 
several  satellite  passes,  between  Nov.1974  and  Feb.  1975. 

This  analysis  should  give  a period  related  to  the  satellite 
spin  and  another  one  related  to  the  Faraday  rotatlon.The  Fa- 
raday effect  shall  not  produce  a imique  period  during  a sa- 
tellite pass,  since  it  depends  on  the  geometry  of  the  magne 
tic  field  and  on  the  TEC;  nevertheless  its  spectrum  is  not 
too  wide  so  that  on  certain  oondltione,  both,  the  Farttday 
and  the  spin  effects  can  be  separated.  Since  the  spin  ef- 
fect depends  on  the  relative  possition  of  transmitter  and 
receiver,  we  have  only  analyzed  nearly  overhead,  N-S  ( mor- 
ning) orblta.  In  fig.  3 the  periods  df  the  analysed  orbits 


Pi^.  Main  periods  of  the  Fourier  analysis  for  different  orbits 
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are  shown  in  a logarithmic  scale*  As  can  be  seen,  in  the 
first  orbits  two  clearly  separated  periods  appear,  the  shor 
ter  one  being  related  to  the  satellite  spin  and  the  longer 
one  to  the  Faraday  rotation*  At  the  beginning,  the  spin  ef 
feet  prevails  and  the  corresponding  period  is  exponentia- 
lly increasing. till  it  reaches  a similar  value  to  the  Fare 
day  rotation  period*  When  this  happens  (about  orbits  496 
or  521)  the  frequency  spectrum  widens  and  the  determination 
of  the  period  of  maximum  amplitude  is  more  difficult*  After 
these  orbits,  the  Faraday  effect  predominates  over  the  sa- 
tellite spin  effect,  and  the  main  period  of  the  Foiirier 
cuialysls  fluctuates  about  a more  or  less  horizontal  line, 
following  the  X£C  variations*  A linear  fit  was  made  till 
orbit  446  to  show  the  exponential  variation. 

The  results  of  this  analysis  indicates  that.  Included 
a generous  security  margin,  the  Faraday  method  to  obtain 
TEC  values  was  applicable,  at  least,  from  Feb.  1975* 

This  analysis  seems  more  accurate  to  determine  the  sa- 
tellite spin  than  the  measure  of  the  temperatures  in  the 
different  components  of  the  satellite*  In  fact,  according 
with  Santullano  (1975)  the  satellite  spin  did  not  produce 
any  anomalous  change  of  themperature  from  December  1st* 

In  the  Faraday  angles  records,  on  the  contrary,  irregular 
rotations  still  appear  after  this  date* 

DETERMINATION  OF  THE  FARADAY  ROTATION* 

To  obtain  the  Faraday  angles  we  used  the  two  close  fre 
quenoles  (40  and  41  MHz)  method*  The  determination  of  the 
first  ooinoidence  of  the  nulls  of  40  and  41  MHz  has  been 
done  by  one  or  two  of  the  three  more  widely  used  methods; 
all  of  them  are  based,  more  or  less  inmediately,  on  the 
well  known  equation 

= 19.75^0 
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that  relates  the  angle  rotated  by  the  polarization  plane  of 
the  wave  of  41  I<iHz(  i'X,.  ) to  the  difference  of  rotation 

Hi 

of  the  polarization  planes  of  the  wave  of  40  lUz  and  41  MHz 

( An  ). 

Method  1. 

This  is  the  most  simple  one.  In  some  records  with  more 
than  10  nulls  of  41  i^Hz,  the  point  where  the  nulls  of  40 
and  41  IslHz  coincides,  can  be  seen  and  determined  with  high 
probability. Besides,  the  null  number  10  from  the  coinciden- 
ce, that  is  recognized  because  of  its  equidistance  from  the 
ad;)acent  nulls  of  40  LiHz,  (cf  .eq.l) , can  also  be  determined. 
If  both  results  are  coherent,  the  determination  is  valid. 

Method  2. 

Prom  eq.l,  values  of  Oai  obtained  for  several 

successive  nulls  with  the  AH  measured  in  the  records.  Prom 
a null  to  the  next  one  the  value  of  0.4]  increases  by  tt 
(or  diminishes  if  the  satellite  goes  to  the  N).  Then, if  we 
substract  Tt  (or  add  if  the  pass  is  S-N)  to  the  value  of 
fX^^vihen  we  pass  from  a null  to  next  one  we  shall  obtain  al 
ways  the  same  value,  the  value  that  corresponds  to  the 
first  analyzed  null.  In  practice  we  often  obtain  different 
values  with  this  procedure  because  of  irregularities  of  the 
ionosphere.  In  these  cases  the  mean  value  is  given  as  the 
true  one  and  the  mean  square  error  as  a measure  of  the  accu 
racy. 

Method  3. 

The  validity  of  eq.  1 relies  on  the  correctness  of  the 
first  order  theory  from  which  it  has  been  deduced.  Howe- 
ver, it  is  foxuid  very  often  in  praetlce  that  the  slope 
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of  the  straight  line  defined  by  eq.  1 in  the  plane  ( , 

) is  different  from  19.75.  To  correct  empirically  for 
the  absence  of  the  higher  order  terms  in  the  theory  , Kers- 
ley  and  Taylor  (1974-)  propose  to  find  the  best-fit  straight 
line,  of  variable  slope,  for  values  of  measured  from 

an  arbitrary  origin  plotted  against  AD.  • Extrapolotion  of 
this  line  to  zero  AD  gives  the  null  of  41  MHz  that  coinci- 
des with  the  null  of  40  MHz  ( A D “ ® ^ • 

In  our  analysis  we  used  Method  1 when  possible; other- 
wise we  used  the  other  two  methods,  that  very  often  gave 
the  same  results.  When  they  did  not,  different  external  cri 
teria  were  used  to  determine  the  most  probable  results. 
Among  these  criteria,  the  rejection  of  negative  or  null  va- 
lues, the  comparison  with  results  of  a near  pass  or  too 

anomalous  variation  of  TEO  within  the  pass,  were  the  most 
useful. 

DATA  QOHSISTMCY. 


An  indication  of  the  accuracy  of  the  data  can  be  obtai 
ned  by  comparing  them  with  the  TEO  values  deduced  from  the 
signals  of  a geostationary  satellite.  An  example  of  such  a 
comparison  with  data  of  the  ATS-P  satellite,  is  shown  in 
fig.  4.  It  corresponds  to  TEC  data  from  the  21st  to  the 
24th  Nov. 1975.  The  T5X)  values  are  plotted  at  the  latitude 
and  local  time  of  their  sublonospheric  points. 

We  drew  the  variation  of  TEC  obtained  from  the  ATS-P 
satellite  on  a plane  at  37SN  and  the  TEC  deduced  from  the 
different  passages  of  the  INTASAT,  in  planes,  perpendicu- 
lar to  this  one, placed  at  the  corresponding  local  time .The 
re  is  a small  variation  of  local  time  from  the  beginning 
to  the  end  of  each  INTASAT  passage  but  we  have  neglected  it 
in  order  to  simplify  the  flgure.The  dashed  lines  in  the  fi 
giure,  correspond  to  extrapolated  values.  As  it  can  be  seen 
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in  the  points  of  the  same  latitude  and  local  time, the  TEC 
values  obtained  from  both  satellites  are  very  similar,  ex 
cept  for  one  oasetthe  pass  at  1930LT  on  the  22nd,aarked 
with  an  arrow  in  the  figure,  where  the  value  deduced  from 
the  Ih'l'ASAT  is  apreolably  higher  than  that  of  the  AIS-F. 
This  anomaly  can  be  caused  by  the  relatively  strong  TEC 
variations  during  this  period,  probably  due  to  a magnetic 
storm  that  started  the  previous  day  with  an  SSC  at  2303nT, 
(corresponding  to  23  30LT)» 

Other  comparisons  were  made  on  several  different  days 
and  similar  agreement  between  the  values  obtained  from 
both  satellites  was  found. 

It  should  be  netted  that  although  these  are  TEC  values 
obtained  at  the  same  local  time  and  at  the  same  latitude, 
they  belong  in  fact  to  different  universal  time  and  are  ob 
talned  at  different  longitudes. 

Another  analysis  was  made  in  order  to  determine  if  the 
apparently  anomalous  variation  of  TEC  for  several  passages 
of  the  INTASAT  was  a real  ionospheric  effect  or,  on  the 
contrary,  was  caused  by  a systematic  error  of  the  calcula- 
tion method.  The  results  seem.;  to  indicate  that  they  are 
in  fact  real  variations  produced  by  irregularities  in  the 
ionosphere. 


DATA  PROCESS. 

As  we  have  already  said  the  INTA:-AT  orbit  was  sun-sln- 
oronouB,  so  that  the  satellite  goes  through  every  latitude 
always  at  the  same  local  time.  Although  the  subionospherio 
point  does  not  follow  the  t>eme  pattern  (as  it  is  obvlons), 
Its  position  is  also  related  to  its  local  time. For  this 
reason,  an  absolute  separation  of  the  latitudinal  and  diur 
nal  effect  on  the  total  electron  content,  is  not  possible 
with  the  kind  of  orbit  of  the  INTASAT.  However,  the  diur- 
nal effect  can  be  diminished  if  a distribution  of  TEC  data 
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Is  made  In  small  Intervals  of  local  time.  We  have  taken  In 
tervals  of  half  an  hour  oentere<i  in  the  hours  and  half 
hours  of  looal  time.  The  data  corresponding  to  each  of  the 
se  intervals f have  been  grouped  by  geographic  latitude  and 
longitude  in  Intervals  of  one  degree  centered  on  the  exact 
degrees*  The  monthly  mean  value  of  all  data  corresponding 
to  the  same  latitude  and  hour  has  been  considered  as  the 
representative  monthly  value  of  that  latitude  and  time* 

The  restrictive  space-time  resolution  that  we  have  in^- 
posed,  causes  the  ntimber  of  TEC  values  at  some  latitudes 
to  be  small,  with  the  result  that  there  can  appear  some 
spurious  peaks  on  the  variation  of  TEC  with  different  para- 
meters, that  do  not  correspond  with  the  true  mean  variation* 
To  avoid  this  difficulty  a small  smoothing  hds  been  perfor 
med  in  latitude  through  a weighted  running  mean  of  three 
values*  In  this  way,  a table  of  TEC  monthly  values  have 
been  obtained  between  30SN  and  442N  for  the  period  Peb* 
1975-Sept. 1976  and  corresponding  to  the  six  different  lo- 
cal time  hours  0900,  0930  and  lOOOLT  in  the  morning  and 
1930,  2000  and  2030  in  the  evening*  Since  several  values 
In  these  tables  were  missing,  mainly  in  the  low  latitudes, 
they  were  cut  to  l»^titudes  between  362  and  432N  for  the 
morning  hours  and  between  362  and  43^  for  the  evening  In 
order  to  obtain  and  homogeneous  space  distribution* 

SEASONAL  VARIATION* 

The  TEC  data  obtained  in  the  way  described  above,  ha- 
ve been  plotted  as  shown  in  figs*  5a-5f«  in  order  to  see 
the  seasonal  variation  for  each  of  the  six  half  an  hour 
periods,  at  different  latitudes* 

As  can  be  seen  from  the  figures,  a olear  consistent 
maximum  in  summer  and  a clear  minimum  in  winter  appear 
In  the  evening  hours*  The  TEC  variation  in  the  morning  is 
slightly  different  for  the  different  hours*  At  lOOOLT  the 
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Fig, 3 a TKC  annual  variation  at  0900  LT 


F ilA!f  J JA  S 0 N DE  F M A MJ  J A S 

1975  1976 


Fig. 5c  TEC  annual  variation  at  1000  LT 


re  Is  a double  oaxlmum  at  the  equinoxes  with  minima  In  win 
ter  (the  deepest  one)  and  in  s\unjaer.  This  variation  that 
appears  more  clearly  at  the  lower  latitudes.  Is  also  appa- 
rent at  0930  LT  although  not  so  well  defined, and  still  less 
distinctly  at  0900LT, 

These  results  can  be  compared  with  the  seasonal  varia- 
tion obtained  at  the  Observatorio  del  Ebro  from  the  BE-B  sa 
tellite  data  (Cardii3,1966,  Gald6n,1970,  Galddn  and  Alberca, 
1971).  In  these  papers  only  midday  values  were  considered. 
The  seasonal  variation  was  similar  to  the  one  found  in  the 
present  paper  for  the  morning  hours  i.e.:two  maxima  at  the 
equinoxes,  a main  minimum  in  winter  and  a secondary  one  in 
summer.  The  amplitude  of  this  variation  increases  with  the 
solar  activity  so  that,  at  the  lowest  region  of  the  solar 
o/ole,  the  equinoxes  maxima  are  only  insinuated.  Since  the 
IIJTASAT  data  correspond  to  a period  of  very  low  solar  acti- 
vity, the  small  amplitude  of  the  equinoctial  rnwrima  is  in 
agreement  with  the  kind  of  variation  found  for  the  BE^B  da- 
ta. 

The  different  variation  of  TEC  found  for  the  different 
hours  can  be  due  to  the  dally  variation,  mixed  up  with  the 
seasonal  variation.  The  morning  hours  correspond  to  the  pe- 
riod of  the  rapid  increase  of  ionization  after  stinrlse.The 
TEO  at  any  time  of  this  period  depends  on  the  time  distance 
to  the  sunrise  moment, and  this  dependence  diminishes  when 
the  time  distance  increase.  From  winter  to  summer,  the  pe- 
riod corresponding  to  the  morning  hours  is  further  and  fur- 
ther away  from  the  siuirlse,  so  that  the  TEC  should  increa- 
se from  winter  to  summer  due  to  this  effect.  The  increase 
is  greater  for  the  earlier  hours.  At  1000  the  daily  varia- 
tion effect  is  not  so  strong  and  the  curves  are  similar  to 
the  curves  found  for  the  midday  period  with  BE-B  satellite 
data. 

Tor  the  evening  hours,  the  daily  variation  also  chan- 
ges from  winter  to  summer.  In  winter,  they  are  far  away 
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from  the  dally  TEC  naxlmumy  that  takes  place  at  midday.  In 
summer  the  maximum  occurs  latter  in  the  afternoon  so  that 
in  the  hourly  period  of  the  data,  the  TEC  is  still  high  at 
this  season.  The  TEC  variation  due  to  this  effect  should 
have  a maxlimim  in  stuomer  and  minimum  in  winter. 


LATITUDIMAL  YARIATIOM. 

In  order  to  obtain  the  latitudinal  variation  we  have 
plotted  in  fig.  6 the  morning  and  the  evening  TEC  values 
between  30S  and  44 -N. 

The  f\ill  dots  correspond  to  all  TEC  values  and  the 
open  circles  only  to  the  homogenised  values  as  described 
above.  We  have  fitted  a straight  line  of  the  form 

TEC  s=.  Up*  a^X  ( A = latitude) 

by  the  least  square  method  to  the  TEC  homogenised  morning 
and  evening  values  of  fig.  6 and  to  the  TEC  values  corres- 
ponding to  each  of  the  six  half  an  hour  periods  mentioned 
before.  The  coefficients  of  the  lines  and  their  most  pro- 
bable errors  are  given  belowt 

LT  ao  e ai  e 

0900  16.51  0.28  - 0.20  0.007 

0930  16.33  0.10  - 0,18  0.003 

1000  16.62  0,23  - 0.17  0.006 

morning  16.45  0.07  - 0.18  0,002 

1930  18.76  0.67  - 0.25  0.017 

2000  18.84  1.33  - 0.27  0.033 

2030  11.48  0.94  - 0.11  0.023 

evening  16.39  0.97  - 0,21  0,024 

As  can  be  seen  the  morning  hours  have  very  similar  values. 
At  the  evening,  the  values  for  2030LT  are  different  from 
the  values  corresponding  to  the  other  periods.  The  errors 
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are  within  the  limits  of  a good  statlscal  significance, The 
difference  between  morning  and  evening  errors  can  be  due 
to  the  fact  that,  depending  on  the  season,  the  ray  path  in 
the  evening  was  some  times  in  darkness,  some  times  partia- 
lly and  some  times  totaly  ilumlneted. 

V/e  can  say  that  the  latitudinal  dependence  of  TEC  Is 
quite  steady,  changing  a little  from  hour  to  hour  In  the 
morning  and  that  it  is  stronger  in  the  evening  than  in  the 
morning  hours. 

NUMERIOAl  MODEL. 


According  with  these  results,  a numerical  mcdel  has 
been  deduce  that  gives  the  TBC  monthly  values  for  the  men- 
tioned six  half  an  hour  periods  as  a function  of  latitude. 

The  model  Is  valid  for  a period  of  low  solar  activity 
similar  to  the  one  analyzed  In  this  report  (annual  mean 
Wolf  niunbers  15.5  for  1975  and  12.6  for  1976)  and  for  a la 
tltudlnal  range  between  36^  and  44SN  for  the  morning  hours 
and  between  38^  and  43-N  In  the  evening. 

A Fourier  time  serie  fit  to  harmonic  4 was  made  separa 
tely  to  each  latitude  and  hour  f6r  the  12  months  Interval, 
so  that,  for  every  hour  and  latitude  we  have 

(TEC);  = €»♦  E C:sin(i^  j ♦f.)  j = U2  (1) 

^ i=1  ^ T i 


where  (TEC)J  means  the  TEC  monthly  mean  value  of  the  j 
month.  This  gives  nine  coefficients  for  every  latitude  and 

hour  ( Cq tA 

For  every  hour,  a third  degree  polynomial  was  fitted 
separately  to  each  Fourier  ooeffiolent  as  a function  of  the 
latltiide  through  the  equation 


(coef),  =K,  ♦ E K,  (lat.-40*r 
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where 


(coef ), 


1=0.4 


lfl-4  " ' = 5.8 

Thus  we  have  a set  of  36  coefficients  fcr  every  hour  that 


will  give  the  mean  value  of  TEC  at  any  given  latitude  (bet 
ween  the  stated  limits)  at  any  month. 


In  table  I we  give  the  values  of  the  Kim  coefficients 


for  the  hours  we  have  TEC  data.  ?/lth  these  values, the  Fou- 
rier coefficients  for  any  given  latitude  can  be  obtained 
through  eq.  2,  for  any  of  the  specified  hours.  With  the 
Fourier  coefficients  the  TEC  for  any  month  of  the  year  can 
be  deduced  through  eq.  1. 


To  give  an  idea  of  the  accuracy  of  the  fittings,  we 
show  in  fig.  7 several  examples  of  the  TEC  annual  varia- 


tion directly  obtained  from  the  data,  and  the  Fourier  fit. 
Me  give  the  curvas  corresponding  to  the  three  morning  and 
evening  hours  for  the  latitude  of  the  Observatory  (41CN). 
As  can  be  seen  the  fit  is  very  good,  although  there  is  a 
small  smoothing  of  several  peaks  due  to  the  fact  that  we 
only  took  the  first  four  harmonics  in  otir  analysis.  Never 
theless,  the  highest  r.  m . s.  TEC  error  that  corresponds 
to  the  1930LT  curve  is  lower  than  0.5.10^5  ^2^ 


The  quality  of  the  fit  of  the  third  degree  polynomials 
can  be  deduced  from  table  II,  where  the  r.m.s.  errors  of 
the  fitted  Fourier  coefficients  are  given.  The  mean  value 
of  the  coefficients  are  also  shown  in  the  table  to  give 


an  idea  of  the  relative  errors.  All  data  are  given  in 
I tuilts  10^^  el  m^.  As  can  be  seen  the  errcrs  are  rather  small 

and  the  higher  one  usually  correspond  to  phases  of  harmonics 
I of  small  amplitude,  so  that  their  influence  in  the  TEC  re- 

f suits  are  still  less  significant. 

I Finally,  in  table  III  we  give  the  differences  between 


the  experimental  TEC  values  and  those  obtained  from  the  mo- 
del. The  monthly  and  latitudinal  as  well  as  the  total,  r.  m, 
s.  TEC  enor  for  every  hour  are  also  given.  As  in  the  pre- 
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56.  8 

-69 

75 

-68 

26 

2 

-28 

49.  0 

45 

49 

43 

33 

19 

27 

33.9 

-69 

79 

-89 

83 

-39 

23 

53.9 

-71 

77 

-67 

35 

-57 

25 

46.0 

-58 

83 

-64 

43 

-46 

20 

43.8 

-36 

56 

-65 

75 

-16 

8 

38.6 

-38 

28 

-46 

52 

4 

-9 

28.  3 

-50 

31 

-17 

-28 

17 

-25 

26.4 

-25 

31 

16 

-43 

21 

-24 

29.  8 

9 

-15 

-4 

13 

8 

-21 

15.  8 

52 

60 

58 

55 

34 

22 

35.  8 

Tabl«  III.-  Differences  between  data  and  model  (x  lO^^el.nT^) 
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Lat.  N 

J 

F 

M 

A 

M 

J 

J 

A 

S 

O 

N 

D 

h 

43* 

54 

-35 

24 

-15 

9 

-10 

22 

-34 

51  -61 

61 

-59 

43.0 

42“ 

25 

-35 

26 

-17 

28 

-45 

42 

-52 

53  -49 

58 

-53 

44.0 

41“ 

39 

-16 

25 

-29 

37 

-49 

55 

-55 

59  -54 

48 

-44 

46.  1 

40“ 

42 

-12 

-4 

-13 

16 

-21 

43 

-41 

52  -65 

46 

-38 

38.9 

39“ 

13 

3 

-34 

37 

-28 

14 

25 

-62 

64  -78 

82 

-47 

49.7 

38“ 

37 

14 

-46 

71 

-62 

25 

14 

-62 

110-122 

112 

-87 

76.  3 

41 

25 

32 

40 

38 

34 

40 

57 

75 

83 

78 

62 

49.4 

20h 


43“ 

23 

-24 

25 

-24 

24 

-20 

16 

-16 

15 

-12 

16 

-18 

20.7 

42“ 

11 

-14 

21 

-23 

22 

-32 

34 

-27 

28 

-34 

23 

-22 

26.4 

41“ 

21 

-16 

9 

-17 

24 

-24 

35 

-44 

40 

-32 

29 

-19 

28.9 

40“ 

29 

-24 

27 

-19 

21 

-19 

20 

-30 

38 

-31 

29 

-22 

27.  5 

39“ 

33 

-38 

36 

-37 

32 

-37 

42 

-45 

45 

-51 

49 

-51 

43.  5 

38“ 

64 

-55 

51 

-49 

55 

-62 

69 

-75 

82 

-80 

76 

-68 

69.  3 

U 

38 

35 

34 

33 

35 

39 

44 

48 

50 

49 

46 

42 

38.  2 

20'30h 


43" 

3 

10 

-17 

17 

-22 

26  -16 

4 

4 

-2 

6 

-8 

14.  2 

42“ 

22 

-13 

-12 

38 

-31 

20  -35 

30 

-8 

-26 

23 

-13 

25.4 

41“ 

20 

-6 

-3 

19 

-42 

49  -33 

30 

-11 

3 

19 

-29 

37.4 

40“^ 

39“ 

27 

-4 

6 

-3 

-39 

62  -31 

4 

-16 

11 

22 

-39 

29.  2 

41 

-38 

-6 

39 

-25 

23  -57 

53 

-6 

-30 

26 

.28 

36.  0 

38“ 

125 

-49- 

110 

84 

26 

31-147 

72 

86 

-58 

-54 

-3 

84.9 

61 

29 

50 

46 

34 

42  75 

44 

40 

32 

31 

26 

41.  2 

(x  lO^^el.in^) 


Tabl«  1II.>  Diffareitees  b«tw««n  data  and  model 


vlous  table  the  units  are  10^^  el  S^.  As  can  be  seen  the  dlf 
ferenoes  and  errors  are  very  small.  Only  7 from  the  504  va- 
lues  reach  10*^  el  m^  and  the  monthly  r.  m.  s.  errors  remain 

~P 

always  below  0,5  lO’*’”  el  We  oan  say  that.  In  general,  the 
errors  of  the  model  are  within  the  errors  of  the  experimen- 
tal data. 

V«ith  TEC  data  obtained  from  the  signal  of  the  BE-B  satel 
lite  recorded  at  the  University  of  Athens,  Klobuchar  (1973) 
deduced  a numerical  model  for  the  mediterranean  area.  The  mo 
del  gives  the  dlxirnal  variation  of  TEC  for  the  foiur  seasons 
of  the  year  and  for  different  solar  activity  levels.  They  go 
from  a solar  flux  In  10,7  cm  of  90  in  spring  to  130  in  win- 
ter (all  In  units  of  ITJ^^  w m^  Hr“^)  the  r.m.s.  errors  for 
the  different  seasons,  that  correspond  to  our  monthly  r,m,s. 
errors,  oscllate  between  0.9  10^^  el  m^  and  1,6  10^  el  m^. 
The  results  of  this  model  cannot  be  properly  compared  with 
ours  because  they  belong  to  periods  of  higher  solar  activity 
and  are  grouped  in  a very  different  way.  Nevertheless,  an 
idea  of  the  oompatlbility  of  the  results  of  both  models  oan 
be  deduced.  The  nearest  conditions  to  our  model  correspond 
to  Klobuchar*s  Season  I (Spring  (march  21) • mean  10.7  om 
flux  equal  to  90  10"^^  w m^  Hz“^),  For  this  period  the  TEC 
values  obtained  from  his  model  between  0900  and  lOOOLf  mre 
slightly  higher  than  our  values, as  expeoted.  Only  at  lOOOLT 
between  36 s and  37BN,  the  values  of  our  model  are  a little 
higher  than  Klobuchar 's.  On  the  evening,  our  values  are 
olearly  higher  than  those  of  Klobuchar.  These  results,  even 
the  last  one  that  seems  a little  suprising  at  first  sight, 
are  in  agreement  with  the  results  gbtalned  by  Oald6n(1968). 
He  finds  that  for  9 and  10  hours,  the  TEC  values  on  March- 
Aprll  1966  are  higher  than  on  Maroh-Aprll  1965  when  the  so- 
lar aetivlhy  was  lower*  At  20  hours  on  the  contrary,  the  va 
lues  of  1966  are  lower  than  those  of  1965*  Also  It  must  be 
noted  that  at  11  hoars,  the  values  of  1965  are  lower  than 
^ose  of  1966.  These  data  correspond  to  41BN.  It  la  possi- 
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ble  that  tha  Tarlatlon  of  TEC  folloirs  the  same  pattern  at 
lower  latitudes,  but  the  inversion  of  differences  in  the 
morning  hours  occur s a little  earlier.  This  would  explain 
the  differences  of  the  two  models  between  36fi  and  37fi  at 
lOOOLT. 

In  fig.  8,  the  isocontours  of  constant  TEC  obtained 
fr<an  the  experimental  data  and  from  the  model  are  shown. 

Is  can  be  seen,  the  model  fit  is  quite  good.  Only  the  re- 
gions of  greater  gradients  of  TEC,  appear  slightly  smoo- 
thed in  the  model.  As  we  said  already,  this  is  due  in  part 
to  the  limitation  of  number  of  Fourier  harmonics  but  also 


to  the  fact  that  a forth  degree  polynomial  would  fit  bet- 
ter than  a third  one  some  of  the  Fourier  ooeffioients.Xhus, 

16  o 

for  instance,  the  isoline  of  10  x 10  el  wr  in  the  0900LT 
map,  that  oovers  the  months  of  March  to  June,  appears  smoo 
thed  in  the  model.  Also,  the  two  tmctima  of  April  and  June, 
as  well  as  the  May  minimum  at  low  latitudes,  do  not  appear 
in  the  model.  However,  all  the  differences  are  small  and, as 
we  said  already,  remain  within  the  error  of  the  experimen- 
tal data. 

Finally,  we  must  point  out  the  limitations  of  the  models 
first  of  all  it  gives  only  mean  values  of  TEC,  is  valid  for 
low  solar  activity  conditions  (mean  Wolf  number  about  14) 
and  fom  the  stated  hours  and  range  of  latitude.  Besides, the 
model  has  been  deduced  for  the  Iberian  peninsula  and  west 
Mediterranean  eons.  The  fitness  to  other  regions  of  diffe- 
rent longitude,  is  limited  by  the  longitudinal  aeyametry  of 
the  ionosphere. 
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